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ABSTRACT 

The origin of dramatically different electron distributions responsible for Comptonization in black hole X- 
ray binaries (BHBs) in their various states is discussed. We solve the coupled kinetic equations for photons 
and electrons without approximations on the relevant cross-sections accounting for Compton scattering, syn- 
chrotron radiation, and Coulomb collisions. In the absence of external soft photons, the electrons are effi- 
ciently thermalized by synchrotron self-absorption and Coulomb scattering even for pure nonthermal electron 
injection. The resulting quasi-thermal synchrotron self-Compton spectra have very stable slopes and electron 
temperatures similar to the hard states of BHBs. The hard spectral slopes observed in the X-rays, the cutoff at 
100 keV, and the MeV tail together require low magnetic fields ruling out the magnetic dissipation mechanism. 
The motion of the accretion disk toward the black hole results in larger Compton cooling and lower equilib- 
rium electron temperature. Our self-consistent simulations show that in this case both electron and photon 
distributions attain a power-law-dominated shape similar to what is observed in the soft state. The electron dis- 
tribution in the Cyg X-1 soft state might require a strong magnetic field, being consistent with the magnetically 
dominated corona. 

Subject headings: accretion, accretion disks - black hole physics - gamma-rays: theory - methods: numerical 
- radiation mechanisms: nonthermal - X-rays: binaries 



1. INTRODUCTION 

The physical processes giving rise to the X-ray/gamma-ray 
emission of accreting black-hole binaries (BHBs) have been 
a matter of debates over the last four decades. The hard-state 
spectra, showing a strong cut-off around 1 00 keV, are well 
described by thermal Comptonization (e.g. lPoutanenlll998l : 
IZdziarski & Gi erlihskil 12004) . while a we ak MeV tail re- 
quires the presence of nonthermal particles (iMcConnell et alJ 
[1994; Ling et aIl[T997h . The origin of seed soft photons for 
Comptonization is, however, much less clear An appar- 
ent correlation between the spectral slope an d the amount 
of Compton reflection (IZdziarski et al.l fl999h argues in fa- 
vor of the a ccretion disk, whil e the observed optical/X-ray 
correlation (iMotch etalJll982t iKanbach et"aL 2001) leans 
toward the synchrotron hypo thesis (e.g. Fabian et al^ 119821 : 
IWardziriski & Zdziarskil200ll) . Interesting questions are then: 
what stabilizes the X-ray spectral slope at a 0.6-0.8, 
and what fixes the temp erature of the emitting plasma at 
fcTe --50-100 keV (ZdziarskietiD II997'; Poutanen' 119981: 
IZdziarski & Gierlifiski 2004)? Do the feedback from the cool 
accretion disk and the thermostatic properties of electron- 
positron pa irs maardt & Mar aschilfT99l iHaardt et al.l[T994l: 
IStern et al.1 1995; Malzac et al.l 1200 lb play a role here? Or 
does the cooling by synchrotron radiation (iNaravan & Yil 
1 1995!) act as a stabilizer? 

In the soft state, BHB spectra are dominated by thermal 
disk emission of temperature kT-QB ~0.4-1.5 keV. At higher 
energies the spectru m is power-law-lik e and shows no signa- 
tures of the cut-off (Grove et al. 1998) extending possibly up 
to 10 MeV (McConne U et al. 2002). This emission is well de- 
scribe d by Comptonization in almost purely nonthermal plas- 
mas (Poutanen & Coppi"1998^■ 'Gierlihski et al."1999t jCoppil 
[l999; Zdziarski et al. 2001; Zdziarski & Gieriihski 200l. We 

' Astronomy Division, Depaitment of Pliysical Sciences, P.O.Box 
3000, 90014 University of Oulu, Finland; juri.poutanen@oulu.li, in- 
drek. vurm@oulu.fi 

- Tartu Observatory, 61602 Toravere, Tartumaa, Estonia; vurm@ut.ee 



can then ask why the electrons are nearly thermal in the hard 
state, and what causes such a dramatic change in the elec- 
tron distribution when tr ansition to the soft state happens. 
iPoutanen & Coppil(ll998l) proposed that the two states are dis- 
tinguished by the way the energy is supplied to the electrons: 
by thermal heating, dominating during the hard state, and by 
nonthermal acceleration, operating in the soft state. However, 
their treatment of Coulomb collisions (using EQPAIR code by 
ICoppil 1 9^2l 1 1 9991) was approximate, and they have neglected 
the effect of synchrotron boiler, involving the emission and 
absorption of synchrotron pho tons, whic h can act as an effi- 
ci ent particle thermalizer (Ghisellini et al.. 19 88). 

iGhisellini et al.l (1 19981) studied for the first time the com- 
bined effect of the synchrotron boiler and Compton cool- 
ing on the electron distribution and photon spectra (but ne- 
glected Coulomb sc attering). They con sidered a two-phase 
corona model (Haardt & Maraschil ll993; Haardt et al. 1994!; 
IStern et al.1 19'95ir where half of the high-energy radiation was 
assumed to be reprocessed by the disk to soft photons. As 
the actual geometry of the emitting region is not known, we 
start from pure synchrotron self-Compton models (i.e. with 
no external soft photons) and compute self-consistently the 
electron (positron) and photon distributions. We then investi- 
gate how the additional soft photons (e.g., associated with the 
inner radius of the cool accretion disk) affect the equilibrium 
distributions and compare the results of simulations with the 
data on Cyg X-1. The preliminary resu lts of this study were 
presented bv lVurm & Poutanen! (l2008ah . 

2. MODEL SETUP 

We consider a black hole of mass 10 Mq, typical for 
stellar-mass BHBs. We assume that the inner accretion flow 
is hot and a lmost spherical , corre s ponding to the advection- 
dominated ('N aravan & Yil 119951: lAbra mowicz et all 119951: 
Esin et al. 199"^ or to the recently dis cover ed luminous 
hot accretion flow solutions ( YuanI 120031 : lYuan & Zdziarskil 
l2004i) . One expects that most of the gravitational energy re- 
lease happens within about R = 10i?s = 3 x 10^ cm (where 
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i?s — 2GM /(? is the Schwarzschild radius) from the black 
hole, and we thus fix the size of the active region at this value 
in most of the simulations. The released energy needs to be 
transferred to electrons via, e.g.. Coulomb collisions with hot 
protons, collective plasma effects, magnetic reconnection, or 
shocks. We assume that the energy transfer to the electrons is 
given by a power-law-injection function dNe/ [dt dj)(x^~'^^'^' 
extending in the Lorentz factor from 7 = 1 to 10'^. To keep the 
Thomson optical depth of the electrons associated with pro- 
tons Tp constant,^ the same number of electrons from the equi- 
librium distribution is removed from the system. In this case, 
the net power is Linj=^i?'^7Vo((7)inj — (7)cq)"ioC^, where 
(7)inj and (7)oq are the mean Lorentz factors of the injec- 
tion function and of the equilibrium distribution, respectively, 
and Nc is uniquely determined by the model parameters and 

(7)cq- 

The injected electrons are cooled by synchrotron emission 
and Compton scattering at timescales much shorter than the 
accretion time. The synchrotron radiation is strongly self- 
absorbed up to hundreds of harmonics, and therefore the cool- 
ing depends strongl y on the high-energy tail of th e electron 
distribution (see e.g. I Wardziriski & Zdziarskill200l]) . The im- 
portance of synchrotron processes is determined by the ratio 
Vb=UbR^c/ Linj, where Ub=B^/ (Stt) is the magnetic energy 
density and ijnj w ^R^cU^sA (so that 77B~3/47rf«0.25 cor- 
responds to an equipartition of the magnetic and radiation en- 
ergy densities, C/B=?^iad)- The seed photons for Compton up- 
scattering can be provided by the synchrotron as well as by the 
external sources, the cool accretion disk being the most natu- 
ral one. The external soft photons are modeled as a blackbody 
of temperature Tbb determined from the Stefan-Boltzmann 
law -Ldisk=47ri?'^o'gBTBB- The cooling by external photons 
depends on the ratio /=idisk/iinj- The total escaping photon 
luminosity is i=Ldisk + ^inj (1 + DUny 

To model self-consistently particle and photon distribu- 
tions, we solve numerically a set of coupled, time-dependent 
kinetic equations for photons, electrons, and positrons de- 
scribing Compton scattering, cyclo-synchrotron emission and 
absorption, electron-electron Coulomb (M0ller) scattering as 
well as pair production and annihilation. Our calculations are 
done in a simple one-zone geometry with a tangled magnetic 
field and isotropic particle and photon distributions. The es- 
cape probability formalism is used to simulate photon escape 
from the region. The detaile d description of the code a nd its 
extensive testing are given in I Vurm & Poutanenl(l2008bl) . 

3. SYNCHROTRON SELF-COMPTON MODELS 

We first assume that the cool disk is sufficiently far away 
and does not supply any seed soft photons to the inner hot 
flow. Thus we consider pure synchrotron self-Compton (SSC) 
models (/=0). We choose the fiducial parameter set L = 
Linj = 10^^ erg s~^ and rp=1.5 (typical for the hard state of 
BHBs. IZdziai-ski et al.|[T997h . i? = 3 x 10^ cm, Tinj = 3 (ad 
hoc), and ?7b=1. The equilibrium electron distribution con- 
sists of a Maxwellian part with fcTc=66 keV and a power-law- 
like tail with the slope modified by cooling To = Finj + 1 = 4 
(where dN^/dj oc j'^"; see the solid curve in Fig.[T^'). The 
synchrotron emission is strongly self-absorbed with only the 
nonthermal tail above 7 >20 contributing to emission above 
the self-absorption energy at >10 eV. As the amount of seed 

' The total optical depth might be larger due to the produced pairs, but for 
parameters considered here, the amount of pairs is negligible. 
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soft (synchrotron) photons is low, the Comptonization spec- 
trum (produced predominantly by the thermal electron pop- 
ulation) is hard with the photon energy index a « 0.9 and 
a cut-off at '^lOO keV, which is similar to the hard state of 
BHBs. A tail produced by single-Compton scattering off the 
power-law electron tail is clearly visible above MeV. 

Variation of the slope of the injected electrons leads to a 
large change in the tail of the electron distribution and dra- 
matic difference in the synchrotron emission. A soft electron 
injection with Finj > 4 leads to efficient thermalization and 
a small amount of soft photons resulting in rather hard ra- 
diation spectra, with the photon energy index a < 0.7 (see 
dashed curves in Fig.[T^). A hard injection gives more power 
to the nonther mal tail and m ore seed photons for Comptoniza- 
tion (see also iGhisellini e t al. 1998; Wardzihski & Zdziarsk^ 
I2OOI . which causes a drop in the electron temperature (see 
the dotted curves in Fig. [TK). A strong 'bump' also devel- 
ops in the tail of the electron distribution at 7 ~3. The syn- 
chrotron emission produced by these electrons is still strongly 
self-absorbed, while the energy losses and ga ins stay close to 
each other for an extended energy interval ( Katarzvnski et alj 
I2OO6I) . In this regime, the ratio of synchrotron heating and 
cooUng rates for a power-law distribution of relativistic elec- 
trons is 7h/7c ~ 5/(rc + 2)."^ Observe that for Fc =3 (i.e. 
for Finj =2) the heating and cooling rates are ba lanced, how- 
ever, such an equilibrium is unstable (lReesll967h . The Comp- 
tonized spectrum for hard injection Finj =2 (see Fig. [1^) is 
much softer than the hard-state spectra of BHBs, even with- 
out any contribution to the cooling from the disk, strongly 
constraining the electron injection mechanism in BHBs. 

The efficiency of synchrotron cooling depends on the mag- 
netic field strength parameterized here via magnetization ryB . 
At small ?7b (see Fig.[T}3'), synchrotron is inefficient and cool- 
ing is dominated by thermal Comptonization. A higher nor- 
malization of the power-law part of the equilibrium electron 
distribution leads to a stronger MeV tail. For 77B ^1 (and 
Tinj > 3), the thermal Comptonization spectrum is very stable 
with a ^0.7-0.9. At large 7/b > 1, the synchrotron thermal- 
ization operates more efficiently and the thermal part of the 
distribution persists to higher energies. The increasing B field 
compensates for the decrease in the power-law tail leading to 
a higher synchrotron emission, which results in softening of 
the Comptonized spectrum. 

Consider now variations of Tp for the fixed Linj . At high 

This expression can be deiived by employing the delta-function approx- 
imation for synchrotron emissivity to calculate the source function and using 
it in t he expression for heatin g by self-absorption (see e.g. iGhisellini etal] 
[T98aiKato"zvnski et aU200^ . 
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Fig. 1 . — Equilibrium photon spectra (upper panels) and electron distributions (lower panels) p^dr/dp (i.e., momentum per log p, where r is the Thomson 
optical depth and p = -^72 — 1 is the dimensionless electron momentum). The results for the fiducial parameter set L = L^^j = 10^^ erg s~^, Finj = 3, 
Tp = 1.5, »7b = 1> / = 0, ai'e shown by solid curves in all panels, (a-a') Dependence on the electron injection slopes: Finj = 2 (dotted curve), 3 (solid) and 
4 (dashed), (b-b') Dependence on magnetization rjg = 0.1 (dotted), 1 (solid), 10 (dashed), (c-c') Dependence on the optical depth Tp = 0.15 (dotted), 1.5 
(solid), 15 (dashed), (d-d') Dependence on the ratio of the external disk photons to the injected power / = 0, 0.1, 0.3, 1, 3 (solid, dotted, dashed, dot-dashed, 
long-dashed curves, respectively) for constant total luminosity L. The electron temperatures and photon spectral indices are given in Table[T] 



Tp, the equilibrium electron temperature drops, leaving fewer 
energetic electrons for synchrotron emission and, therefore, 
reducing the number of seed photons for Comptonization 
(Fig. [T];). This, in turn, results in the harder photon spectra 
produced by saturated Comptonization (by thermal electrons) 
and a weak high-energy tail (produced by nonthermal elec- 
trons), ver y similar to the ultrasoft spectra of BHBs (see figs. 
8 and 9 in lZdziarski & Gierliriskiir2004) . At smaller Tp, the 
higher electron temperature leads to a stronger synchrotron 
cooling and to a lower Comptonized luminosity and, there- 
fore, softer Comptonized spectra. 

Let us now apply the developed model to the hard state 
of Cyg X-1. The Me V tail observed there with aucV ~ 2 
dMcConnell et al J 120021) constrains the injection slope to be 
Finj < 2aMcV = 4. Then the hard X-ray spectra with a <0.7 
and a high-energy cutoff at ~100 keV (see Fig. |2]) require 
Tp ^ 1 and low r]B 5,0-l (se e also IWardziriski & ZdziarskH 
120011; iMcConnell et al.ll2002l) . Any additional soft photons 
from the disk will make the spectrum softer, reducing 7/b even 
more. The low magnetic field rules out magnetic reconnec- 
tion as the energy dissipation mechanism. This also implies 
that electrons cannot be thermalized by the synchrotron self- 
absorption. On the other hand, if the size of the active region 
is i?~60i?s. Coulomb scattering becomes important (as its 
influence grows linearly with size for constant L, see, e.g., 
ICoppi 1999; Svensson 1999), and it can thermalize electrons 
at a rather high tem perature of fcTe^lOO keV a s observed in 
Cyg X-1 ( Gierlinskieral.lll997tlPoutanenl[T99l 

We reiterate that the whole spectrum here is produced by 
the SSC mechanism.^ Its thermostatic properties fix the elec- 
tron temperature at 50-100 keV (for Tp ~ 1) and stabi- 

^ ADAF-based model s also consider SSC as the main cooling mechanism; 
see INaravan et all fT998h for the review. 



lize the spectral slope at a^O.7-0.9. The feedback from the 
disk (Haardt & Maraschi 1993; Haa rdt et al.ll994tlStern et all 
il995i;.Malzac et al..,2001.) does not seem to be needed. 

4. SPECTRAL TRANSITIONS AND THE ROLE OF DISK 
PHOTONS 

The spectral transitions observed in BHBs are most prob- 
ably accompanied by a change in the geometry of the ac- 
cretion disk. The cool outer disk moves toward the central 
black hole causing an increasing flux of the soft photons to 
the central h ot flow (Esinetal. 1997; Poutanen et al. 19971 
iPoutanen & Coppill 19981) . which we simulate here by increas- 
ing / (see Fig. [TJl). Higher soft photon flux leads to faster 
Compton cooling and lower equilibrium electron temperature, 
making the nonthermal part more pronounced. The resulting 
photon distribution changes from the hard, thermal Comp- 
tonization dominated, spectrum to the one dominated by the 
disk blackbody, with a nonthermal tail extending to tens of 
MeV, which becomes harder at higher /. The spectral changes 
triggered by varying / are similar to the one observed in Cyg 
X-1 (see Fig. |2]l. A detailed comparison with Cyg X-1 spectra 
shows, however, that other parameters change too. 

Compared to the hard state, the soft state corresponds to a 
higher total luminosity. The MeV tail is harder aMeV ~ 1-6 
(iMcConneU et al.ll200 2'). and therefore Finj < 3.2. If the tail 
of the black-body at 3-10 keV (see Fig. |2]i is produced in the 
same emission region, it requires a rather hot thermal popu- 
lation of electrons, which needs high tjb for the synchrotron 
thermalization to operate (because Coulomb thermalization is 
not efficient under the conditions of strong Compton cool- 
ing). This would be consistent with the magnetically domi- 
nated emission region. Alternatively, there may be additional 
heating mechanisms operating. Also the tail might be a re- 
sult of Comptonization in the hot ionized skin of the disk, not 
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Fig. 2. — Spectral st ates of Cyg X -1. Crosses show the unfolded 

spectral data presented by IZdziarski et al.l j2002l) . The model spectra with 
interstellar absorption (described by the column density A^h) and Comp- 
ton reflection (described by the solid angle Q and ionization parameter ^, 
IMaedziarz & Z dziarski 1995) taken into account are shown by the solid 
curves. The parameters for the hard-state model are: L = 2.7 X 10"^^ erg 
s-i, R = 1.8 X 10** cm, / = 0, T,^^ = 3.8, Tp = 2.5, rje = 0.083, 
ATjj = 3 X 10^1 cm-2, n/(27r) = 0.2, ? = 0. The soft state can be 
described by L = 4.85 X lO^'' erg s"!, R = 3.85 X lO'' cm, / = 2.13, 
Tinj = 2.2, Tp = 0.3, ??B = 1.5, ATh = 5 X lO^i cm'^, f2/(27r) = 0.7, 
i = 100. 

directly related to the emission we discuss here, but this in- 
terpretation might not be easily reconciled with the fact tha t 
the disk is stable, while the tail varies (Chur azov et al.ll2"001|) . 
While the dramatic changes in the electron and photon dis- 
tributions between the states are mainly caused by variations 
of the disk luminosity, it is obvious that other parameters do 
change during the transition. We stress that none of the pre- 
sented models requires any ad ditional thermal hea t ing, w hich 
is different from the in odels of iPoutanen & Coppil (Il998l) and 
lGierhnskietan(ll999h . 

5. CONCLUSIONS 

The hard state of BHBs can well be described by the quasi- 
thermal SSC mechanism. The feedback from the cool disk 



is not needed to stabilize the spectral slope and the electron 
temperature. Electrons can be injected to the active region 
with the power-law spectrum, but Coulomb scattering and 
synchrotron self-absorption thermalize them efficiently. This 
reduces the need for mysterious 'thermal heating' that was in- 
voked previously to explain thermal Comptonization spectra 
of BHBs. The MeV tail together with the hard X-ray spectra 
of BHBs with photon indices a <0.7 and a cutoff at 100 keV 
require rather low magnetization r/B < 0.1 and a large size of 
R > 60i?s- In that case, magnetic reconnection can be ruled 
out as a source of energy. 

At high optical depth of the emitting region Tp > 10, in the 
absence of disk radiation, the spectrum is close to saturated 
Comptonization, peaking at a few keV. This Wien-type spec- 
trum might be associated with the ultrasoft state of BHBs. At 
low Tp, the electrons are hotter and the spectra are softer due 
to the efficient synchrotron cooling. 

A behavior similar to what is observed during the spectral 
state transitions in BHBs can be reproduced by varying the ra- 
tio of injected soft luminosity and the power dissipated in the 
hot phase, which could be caused by varying the radius of the 
inner cool disk. The increasing Compton cooling causes dra- 
matic changes in the electron distribution from almost purely 
thermal to nearly nonthermal. The photon distribution also 
changes from quasi-thermal SSC to the nonthermal Comp- 
tonization of the disk photons. In the soft state of Cyg X- 
1, a strong magnetic field can thermalize electrons at suffi- 
ciently high temperature, which is consistent with a magneti- 
cally dominated corona being responsible for the high-energy 
emission. 
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